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Abstract 

We present an exact analytical solution of the Einstein equations 
with cosmological constant in a spatially flat Robertson- Walker met- 
ric. This is interpreted as an isotropic Lemaitre-type version of the 
cosmological Friedmann model. Implications in the recent discovered 
cosmic acceleration of the universe and in the theory of an inflationary 
model of the universe are in view. Some properties of this solution are 
pointed out as a result of numerical investigations of the model. 
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1 Introduction 



Recent astrophysical investigations 0, || demonstrate that the expansion of 
the universe is accelerating rather than slowing down. This may change our 
picture of the universe suggesting that we live in an accelerating and flat uni- 
verse. However, from the theoretical point of view intensive efforts are done 
M in order to accommodate the relativistic cosmology with an accelerating 
universe. One of the major results of this situation is the fact that the most 
popular solution (inspired by what is happening in the theory of inflationary 
cosmology) is to consider models which satisfy Einstein equations with cos- 
mological constant in order to induce cosmic acceleration. Actually, several 
models until now not being in the main stream of the modern cosmology, are 
again in view of cosmologists. We believe that it is necessary to more care- 
fully investigate new possible versions of the the Friedman-Robertson- Walker 
(FRW) model J|, |6j in order to see how may we implement here cosmic ac- 
celeration without loosing the well-known results of the model in describing 
early stages of the universe. 

Another important problem is to try to find simple analytically solvable 
models which may be easily interpreted from the physical point of view. 
Moreover, it is clear that many further developments of the more complicated 
models have to be done using algebraic or numerical methods on computers. 
Then it is important to have at least one simple analytically solvable model 
for testing the computational methods and verify their degree of confidence. 

This article is dedicated to the revealing of an analytical solution of the 
Einstein equations with cosmological constant in a version of the FRW model 
with flat space metric ||. We present some properties of this solutions and 
the possible implications to the more accurate description of the inflationary 
and accelerated universe. 

We must specify that the model with cosmological constant and k = 1 
is known in literature as Lemaitre-FRW model but, unfortunately, this does 
not have analytical solution. In other respects, we believe that the solution 
for k = we present here is ignored by the investigators of this cosmology 
since they are forced to find solutions of the Einstein equations in terms of 
the measurable quantities (as the Hubble constant/parameter, the redshift 
or the deceleration parameter) and searching for an approximate solution 
for the scale factor of the universe as a Taylor expansion in terms of these 
quantities. Thus the existence of an analytical solution is hidden by the large 
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number of solutions available only for restricted periods from the history of 
universe. 



2 The model and its solution 

Let us consider the Lemaitre-type version of the FRW model with flat space 
(i.e., k — 0) 0, positively defined cosmological constant, A, and linear depen- 
dence between the density of energy, e, and the pressure, p. In general, in an 
arbitrary local chart (or natural frame) of coordinates x^ (fi, u, .. = 0, 1, 2, 3), 
the Einstein equations of this model (written in natural units with c = 1), 

r; - - ^; = ^gt; , (i) 

involve the usual classical stress-energy tensor 

t; = ( e + P )«x - P s; , (2) 

which depends on the covariant four-velocity = dx^/ds. 

In the preferred frame with Cartesian coordinates x° = t, x which has 
the usual FRW line element with k = 0, 

ds 2 = dt 2 - a(t) 2 dx 2 , (3) 

we have u° = u = 1 and u % = (i = 1, 2, 3) such that the Einstein equations 
reduce to the simple system 

/An 2 

= A + 87rGe (4) 
= A-8nGp (5) 

where we denoted ' = d/dt. Assuming that e and p are functions only on t 
and satisfy a linear equation of state, namely 

p = ne, (6) 

one can integrate this system with the conditions 

e(0) = e , a(0) = l, (7) 
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at the present time t = 0. Indeed, 
equation 

d / a\ k + 1 
(it Va, 



(|) and (|6|) are equivalent with the 
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A-3 



(8) 



which, after a little calculation and exploiting (0), leads to the final solution 
a(t) = cosh — - — [K + l)t 

I / 1 2/3(k+1) 

(9) 
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while the density of energy reads 

e(t) = e a(t)- 3{K+1) . 



(10) 



Finally, we calculate the Bang time, —to, from the initial condition a(— to) = 
as 

2 , /" 87rGe 
= — — argcothW H — , (11) 



3A(k + 1 



completing thus the solution of this simple model. 

Let us observe that the conservation law = is implicitly fulfilled 
since ( PH| ) is simultaneously the solution of this conservation equation which 
takes the form 

(12) 



e + 3e(K + l)- = 
a 



when the equation of state (^) is accomplished. Moreover, it is important 
that the equation (|T0|) is independent on A since then it holds even in the 
particular case of A — >• when we recover the usual solution of the FRW 
model with k = §, |, 



a(t) 



1 + J6irGeo(K+ l)t 



1 2/3(k+1) 



(13) 



The interesting new features of our model are due to the cosmological 
constant which gives an hyperbolic character to the general model or to its 
asymptotic behavior. It is clear that the universe devoid of matter becomes 
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a de Sitter one but it is remarkable that in the far future the time behavior 
of a(t) is also of the de Sitter type since for very large t we have 



a(t) 



~ a e 



a 



8nGe 
A 



2/3(k+1) 



(14) 



In other words if we rescale the Cartesian space coordinates with the factor 
a then the line element (BI) becomes just the de Sitter one in the limit of 
t — > oo. Thus we can conclude that our solution is compatible with the 
well-known models describing the Big-Bang scenario and for late behavior it 
shows off characteristics of a de Sitter spacetime. 

Many physical effects depend on the form of the Hubble function that in 
our model reads 



H(t) 



1 



3(A + 87rGe ) 



A 



87rGe 



1 + V 1 + ^f a tanh^(/t + 1) t 
Obviously, in the late time limit we have 

lim Hit) = 

while the actual value of the Hubble function (i.e., the Hubble constant) 



• (15) 



#o = Hit) 



I A + 87rGe 



t=o- 



(16) 



depends on the actual value of the density of energy, e , and A. 



3 Numerical investigations 

This section is dedicated to several numerical results we obtained using our 
analytical solution which reveals some interesting features of the model. First 
we have plotted the behavior of the scale factor in time for a constant value 
of the cosmological constant A. In this purpose we considered an accepted 
value of the actual mass-density of the universe as po — 3 ■ 10~ 28 Kg/m 3 
(in IS units where e = Po c 2 ) an d of the actual Hubble constant as H = 
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Figure 1: Evolution of the scale factor a(t) from the Big-Bang (left panel) 
and the fast growing of the scale factor at late time t = 10 18 ...5 10 19 s (right 
panel) 

(1.8 10 10 years)~ L [H]. Then, using (|TB|) we can estimate the cosmological 
constant as being A m 0.977 10 -52 m~ 2 which will be the value used in the 
next numerical investigations. We also choose k = —1/2. Two first graphs 
are showed in Figure [I] for different intervals of time, one before the actual 
time t = and the other one for late behavior of the universe. It is obvious, 
especially from the second picture (right panel) the fast growing of the scale 
factor for late time, showing the accelerated stage of the universe described 
by this model. We have to mention that we estimated the Big-Bang time to 
using the above numerical values and formula (jTTf) . Thus we have used in our 
figures the resulting value t « 0.1665 10 19 s (~ 52.7 billion years). This value 
is larger than the recognized one and appears here due to the small value 
of the mass-density we used. A greater value for p, including dark matter 
and other unknown components, (for example one ten times larger) gives a 
Big-Bang time to ~ 29 billion years. . 

Figure [| presents the time behavior of the Hubble constant (the left panel) 
as defined in fll5|) for the same parameters as for the above graphs. As we 
pointed out in the comments ending the last section, the Hubble function 



is evolving to a constant value yA/3 specific for the de Sitter model. This 
conclusion is coupled to the late time behavior of the scale factor (see above) 
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Figure 2: Evolution of the Hubble function (left panel) and of the redshift 
defined as l/a(t) (right panel) 



which is also specific to a de Sitter metric. In the right panel of the same 
Figure |2] we plotted the redshift defined as l/a(t). 

Next plots, Figure [| are dedicated to the time behavior of the acceleration 
defined here as the second time derivative of the scale factor, a(t). Before the 
actual time t = (presented in the left panel) shows the decreasing of the 
acceleration after the Big-Bang at t > —to- In the right panel we observe the 
increasing of the acceleration at later time in the universe evolution. Thus 
we can conclude that it will be possible to use our solution for modeling the 
so-called "cosmic-acceleration" in discussion in the modern astrophysics. 

Last of this type of plots, Figure |4], represent the time behavior of the 
density function e(i) given by equation (PH|) and of the deceleration defined 
^■a{t)/\a{t)H{tf} HEg. 
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Figure 3: Evolution of the acceleration, a(i), before the present time (left 
panel) and after t = (right panel) 
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Figure 4: Evolution of the density before the actual time t — (left panel) 
and of the deceleration, a(t)/[a(t) H(t) 2 ], after the actual time (right panel) 
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